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Static-light mesons on a dynamical anisotropic lattice 
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Abstract. We present results for the spectrum of static- light mesons from Nf = 2 lattice QCD. These results 
were obtained using all-to-all light quark propagators on an anisotropic lattice, yielding an improved signal 
resolution when compared to more conventional lattice techniques. In particular, we consider the inversion 
of orbitally-excited multiplets with respect to the 'standard ordering', which has been predicted by some 
quark models. 
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1 Introduction 

It has long been known that hadrons containing a single 
charm or bottom quark exhibit approximate heavy quark 
spin and flavour symmetries and can be described to a 
reasonable approximation by an idealised system in which 
the heavy quark is taken to be infinitely massive pQ. In 
this limit the heavy quark symmetries become exact and 
the system consists of light degrees of freedom bound to a 
static colour point source. The static approximation can 
be systematically improved using a power-counting rule 
in ylQCD/niQ, which forms the basis for heavy quark ef- 
fective theory (HQET). For hadrons containing a bottom 
quark this approach is particularly suitable. In this case 
one expects leading-order corrections to be already quite 
small, at the level of ten percent. 

On the lattice, mass-dependent errors which arise from 
the discretisation of the Dirac operator can limit the ac- 
curacy of simulations of heavy-light systems. In principle, 
one can avoid such errors by using the HQET lagrangian 
to describe the heavy quark dynamics. An alternative ap- 
proach is to perform a number of simulations both at 
unphysically light heavy quark masses using the lattice 
Dirac operator and in the static limit and to interpolate 
the results to the physical heavy quark mass guided by 
the predictions of HQET. Therefore, precise calculations 
in the static limit are of considerable interest to the lat- 
tice community. In addition, little is known experimen- 
tally of the excited state spectrum of heavy-light hadrons 
and accurate simulations of their simpler static-light coun- 
terparts will have considerable phenomenological impact. 
However, traditionally, simulations of static-light systems 
have been extremely noisy. Ultimately, this stems from the 
fact that the use of conventional point-to-all propagators 
for the light quarks restricts interpolating operators for 
the static-light hadron to a single spatial lattice site. 



This problem could be overcome if it were possible to 
compute all elements of the light quark propagator. This 
would allow source and sink operators for the static-light 
correlators to be placed at every spatial lattice site, yield- 
ing a dramatic increase in statistics. Indeed, simulations 
of static-light mesons are the natural testing ground for 
all-to-all techniques since they require a single light quark 
propagator per configuration while the static quark prop- 
agator is easy to evaluate. It is therefore not surprising 
that a number of groups have attempted to tackle this 
particular problem [21EJ - 

In this contribution, we describe a preliminary study of 
the excited state spectrum of static-light mesons using the 
all-to-all propagator technique introduced in Ref. 0]. In 
this approach, the important contribution of the low-lying 
eigenmodes of the Dirac operator to the fermion propaga- 
tor is computed exactly. A noisy estimate is then used to 
correct for the contribution of the remaining eigenmodes. 
This estimate is obtained by inverting the Dirac operator 
on a set of Z4 noise sources which have been partitioned 
in some subset of quarkfield indices. This partitioning or 
'dilution' yields a substantial reduction in the variance of 
the stochastic estimate. This study has been performed 
on a 3+1 anisotropic lattice, where the temporal lattice 
spacing a t is much finer than the spatial lattice spacing 
a s . The fine temporal lattice spacing proves particularly 
useful when determining a plateau in the effective masses 
of the higher excited states. The correlators have been 
computed on Nf = 2 background configurations. 



2 Interpolating operators 

Due to heavy quark spin symmetry, states which lie in the 
same hyperfine multiplet are degenerate and in the con- 
tinuum it is conventional to label these multiplets by the 
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we tested a number of operators across the lattice repre- 
sentations. Further details of our choice of interpolating 
operators will be presented in Ref. 0. 



Table 1. Irreducible representations of the group of 24 proper 
spatial lattice rotations with low-lying angular momentum con- 
tent. 

angular momentum and parity of the light degrees of free- 
dom pj. For example, in the static limit the pseudoscalar 
and vector correspond to Je Pe = \ ■ 

At finite lattice spacing, it makes sense to classify states 
according to the irreducible representations (irreps) of the 
octahedral point group Oh [E]- This 48 element group 
is the direct product of the 24 proper spatial rotations 
allowed by the lattice with the group consisting of the 
identity and the space inversion operator. The relation- 
ship between these representations and continuum quan- 
tum numbers is not difficult to determine. Restricting the 
continuum irreps to the elements of the lattice symmetry 
group generates representations for Oh which are in gen- 
eral reducible and identification of the constituent irreps is 
a straight-forward exercise in group theory. Here, we take 
advantage of the degeneracies which arise from the heavy 
quark spin symmetry and construct interpolating opera- 
tors specifically for the light degrees of freedom. Therefore, 
our operators transform according to the double- valued or 
spinorial representations of Oh- There are 6 such represen- 
tations 

Giu, Gig, G2u, G2 g , H u , H g . (1) 

where the subscripts u and g denote representations which 
are odd and even under spatial inversion respectively. 

Tabled lists the irreps of the subgroup of proper rota- 
tions and the angular momenta of their lower-lying con- 
stituent states. From this table we see that states with 
the same angular momentum but with different J z may 
be scattered across the lattice irreps. For example, states 
which transform according to the rows of the 6-dimensional 
spin | representation are divided between the G2 and H 
irreps. In a numerical study we should therefore expect to 
observe energy levels in these irreps which are degenerate, 
up to lattice artifacts, corresponding to the Ji = \ states. 
In this study, we were particularly interested in orbital ex- 
citations of the static-light meson. There is a single S-wave 
state with J^ e quantum numbers \ , the P-wave doublet 
is i + , | + and the D-wave states are labelled | , | .Of 

these, only the S-wave and the i + P-wave which lie in the 
Gi u and Gi g irreps can be accessed using local operators, 
while all other irreps require the use of extended interpo- 
lating operators. It is worth noting that when one uses 
all-to-all light quark propagators the construction of ex- 
tended interpolating fields requires no additional fermion 
matrix inversions and the extra computational cost in- 
curred is minimal. In order to determine accurately the 
excited state spectrum it is important to choose opera- 
tors which couple strongly to the states of interest and 



3 Simulation details 

The lattice Dirac and gauge actions used in this study have 
previously been described in detail in Refs. Both 
are specifically formulated for 3+1 anisotropic lattices. 
The Wilson-like fermion action is accurate to C(at,af) 
and employs stout-smeared spatial links to reduce radia- 
tive corrections. The gauge action incorporates tree-level 
Symanzik and tadpole improvement and is accurate to 

The anisotropic lattice breaks hypercubic symmetry 
and introduces two new parameters: the bare quark and 
gauge anisotropies, denoted £° and which must be 
tuned such that the measured anisotropy, a s /at, assumes 
a consistent value. In dynamical QCD the tuning proce- 
dure is quite complicated and is detailed in Ref. In 
this study we use bare anisotropy values taken from that 
paper of £° = 7.52 and £° = 8.06 which yields a measured 
anisotropy of 6 with a three percent error. 

We employ the standard Eichten-Hill action for 
the static quark and the corresponding static quark prop- 
agator is a product of unsmeared temporal links. 

Measurements have been performed on 249 background 
configurations on an 8 3 x 80 lattice with a spatial lattice 
spacing of about 0.2fm. The mass of the sea quarks and 
the light valence quark was found to be approximately the 
strange quark mass. For the light quark propagator, we 
computed the contribution of 100 low- lying eigenvectors 
of the Dirac operator. Two independent Z4 noise sources 
which were diluted in time and colour were used in the 
stochastic estimate. 

To extract excited state energies we applied the stan- 
dard variational analysis to correlation matrices which 
were obtained by applying 5 or 6 levels of Jacobi smearing 
to the light quark fields. This was facilitated by the use 
of all-to-all propagators which meant that no additional 
inversions were required to smear the quark fields. 



4 Results 

Fig. n shows effective masses for the lowest-lying states 
across a number of the lattice irreps. In this case the lowest 
effective mass corresponds to the S-wave \ state and the 
higher- lying states are orbital excitations. The highest- 
lying states lie in D-wave channels. Clear plateaux are 
evident in each channel. 

In addition, we obtain strong signals for a number 
states in each of the lattice irreps. Fig.|21shows the ground 
state and the first two excited states in the Gi u (S-wave) 
irrep. These excited states shown here are almost certainly 
radial excitations in the S-wave channel. 

Since we are working in the static limit the energies 
determined by fitting to the effective masses contain an 
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Fig. 1. Effective masses for the lowest-lying states in 5 lattice Fi S- 3 " Mass differences between the ground state S-wave mul- 



irreps. These include P-wave and D-wave excitations. 
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Fig. 2. Effective masses for the lowest-lying states in the Gi u 
(S-wave) irrep. 



5 Conclusions and outlook 

We have presented preliminary results for the spectrum of 
static- light mesons in Nf = 2 QCD. We have been able to 
compute energy differences between the S-wave multiplet 
and a number of orbital and radial excitations. The use of 
all-to-all propagators allows us to minimise the statistical 
uncertainty in our measurements and facilitates the use of 
spatially-extended interpolating operators and the varia- 
tional approach. The results presented here have allowed 
us to make some tentative remarks about the nature of 
the spectrum; however, we have not yet made a serious at- 
tempt to identify the continuum quantum numbers of the 
observed states. More recently we have performed similar 
simulations using larger operator bases on two different 
spatial volumes. These runs should allow us to assess fi- 
nite volume effects and identify multi-particle states. We 
are currently analysing the results of these simulations and 
we intend to present our findings in the near future F7|. 



unphysical shift. This shift cancels in energy differences, 
which are therefore physically meaningful. Fig. [3] plots 
the energy differences between the lowest-lying state in 
the Gi u representation (i.e. the S-wave multiplet) and a 
number of excited states. We have not presented these 
preliminary results in physical units because we have per- 
formed our simulation on quite a small spatial volume 
which might distort the energies of the orbitally and ra- 
dially excited states. The main point here is the precision 
to which we can compute the spectrum. Nevertheless, it is 
interesting to note the position of the lowest-lying states 
in the H u and G2 U irreps. Naively, one expects the light- 
est state in the H u representation to correspond to the 
| D-wave while the | ought to be the lightest state in 
the irrep. However, there appears to be no significant 
splitting between these states which may indicate that the 
| multiplet is lighter than the § . Such a scenario has 
been predicted by quark models |12I13| . 
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